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ABSTRACT 

The Solar Diameter Monitor measured the duration of solar meridian transits during 
the 6 years 1981 to 1987, spanning the declining half of solar cycle 21. We have combined 
these photoelectric measurements with models of the solar limb-darkening function, de- 
riving a mean value for the solar near-equatorial radius of 695.508 ± .026 Mm. Annual 
averages of the radius are identical within the measurement error of ± .037 Mm. 

Subject headings: Sun: diameter, helioseismology 
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1. Introduction 

The Sun is the only star for which reasonably pre- 
cise values of the mass, surface radius and luminosity 
are known. The solar mass Mq is known from plan- 
etary motion, with accuracy limited only by the un- 
certainty in the gravitational constant G. The solar 
radius can in principle be obtained from direct opti- 
cal measurement of the solar angular diameter, given 
the very accurate determinations of the mean distance 
between the Earth and the Sun. In solar modeling, 
the value R Q = 695.99 Mm (Allen 1973) has been 
commonly used. The models are calibrated to this 
photospheric radius, in the present paper defined by 
the point in the atmosphere where the temperature 
equals the effective temperature, by adjusting some 
measure of the convective efficacy, such as the mixing 
length. 

Recent accurate observations of solar f-mode fre- 
quencies from the SOI/MDI instrument on the SOHO 
satellite (e.g. Kosovichev et al. 1997) have raised 
some doubts over this value of Rq. The frequen- 
cies of these modes are predominantly determined 
by GMq/Rq. By comparing the observed frequen- 
cies with frequencies of solar models calibrated to 
R Q = 695.99 Mm Schou et al. (1997) and Antia 
(1998) concluded that the actual solar radius was 
smaller by about 0.3 Mm than the assumed radius 
of the model. Other aspects of the modeling of the 
solar f modes may affect their frequencies at this level 
(e.g. Campbell & Roberts 1989; Murawski & Roberts 
1993; Ghosh, Antia & Chitrc 1995). Thus it is obvi- 
ously important to obtain independent verification of 
the proposed correction to the solar radius. 

There are indeed significant uncertainties associ- 
ated with the currently adopted radius value. These 
are related to the problem of the definition of the so- 
lar limb adopted in the radius determinations, and 
the reduction of the measured value to the photo- 
sphere. It is not clear how the value quoted by Allen 
(1973) was obtained. However, it appears that the 
more recent determinations, which are generally con- 
sistent with Allen, in most cases refer to the inflection 
point of the solar limb intensity. According to solar 
atmospheric models this corresponds to a height of 
about 0.3 Mm above the photosphere, thus perhaps 
accounting for the radius correction inferred from the 
f-mode frequencies. 

The uncertainty in the precise definition of the 
measured values of the solar radius highlights the need 



to combine the observations with careful modeling of 
the quantity that is observed. Here we consider a long 
series of observations obtained with the High Altitude 
Observatory's Solar Diameter Monitor (Brown et al. 
1982). This is based on a definition of the solar limb 
which minimizes the effect of seeing (Hill, Stebbins & 
Olcson 1975). By combining daily data obtained over 
more than 6 years, extending between solar maximum 
and solar minimum, the possible effects of solar activ- 
ity can be checked. The analysis of the data is carried 
out by means of a model of the solar limb intensity, 
following as closely as possible the actual procedure 
used in the reduction of the data and testing for the ef- 
fects of seeing. In this way we have eliminated several 
of the uncertainties affecting earlier determinations to 
arrive at what we believe to be an accurate measure 
of the solar photospheric radius. 

2. Observations 

The Diameter Monitor instrument and its associ- 
ated observing procedures were described in detail by 
Brown et al. (1982); it operated between August 1981 
and December 1987. It consisted of a meridian-transit 
telescope arranged to allow the solar image to drift 
across a fixed detector package each day at local noon. 
A filter system confined the bandpass of the observed 
light to a 10 nm band near 800 nm. The horizontal 
(east/west in the sky) diameter was obtained by tim- 
ing the passage of the solar limbs across each of two 
linear detector arrays aligned end-to-end. Each detec- 
tor pixel subtended 1 1T in the sky in the direction of 
the apparent solar motion and 80 " in the perpendic- 
ular direction. In addition, the vertical (north/south 
in the sky) diameter was measured, although less pre- 
cisely. For the purposes of this paper we shall there- 
fore consider only the horizontal diameter. 

An automated guiding system assured that the so- 
lar disk transited the detector arrays centrally, so that 
a true diameter was measured. Each readout of a 
detector array yielded a sample of the solar limb- 
darkening function; by reading the detectors at a 
32 Hz rate, the instrument obtained samples at inter- 
vals comparable to the seeing-change time. The in- 
strument applied a real-time edge-finding algorithm, 
and stored the resulting edge positions. This process 
was performed for the transits of both west and east 
solar limbs, so a transit duration could be measured. 
Ancillary quantities were also measured each day, in- 
cluding seeing and scattered-light parameters. 



2 



The edge-finding algorithm used was the Finite 
Fourier Transform Definition (FFTD) described by 
Hill et al. (1975). The procedure involves forming 
the convolution of the observed limb-darkening func- 
tion with a set of weights that are nonzero only within 
a certain window of width a. The edge was then de- 
fined to be the position of the center of the window 
for which the convolution crossed zero. 

The FFTD has two important features. First, by 
a suitable choice of weights one can eliminate the 
first-order sensitivity of edge position to seeing, for 
some chosen width of the seeing point-spread func- 
tion. Daytime seeing is commonly both poor and 
variable, and inflection-point definitions of the limb 
position are highly sensitive to this variability. Sec- 
ond, the FFTD depends on a free parameter, namely 
the window width a. Applying the FFTD to the solar 
limb-darking function yields diameters D(a) approx- 
imately equal to 

D(a) = D - aa , (1) 

where D is the true angular distance between the 
nearly discontinuous intensity jumps at east and west 
limbs, and a is proportional to the intensity gradient 
in the last 25 B inside the limb, and hence to the ver- 
tical temperature gradient in the photosphere above 
an optical depth of roughly 0.2. By measuring the 
diameter over a range of a and extrapolating the re- 
sults to zero window width, one can obtain a value 
that is largely independent of changes in the slope of 
the limb-darkening function. Since this slope proves 
to vary significantly during the solar cycle, using mul- 
tiple window widths is necessary to measure the slope 
and remove its effect from the raw diameter measure- 
ments. 

Data were taken with the Diameter Monitor on 
any day for which a successful observation seemed 
possible; many of the observations were therefore 
corrupted by clouds or (less often) instrument fail- 
ures. Accordingly, we used several selection criteria 
to choose the observations to be used for the analysis, 
in the end retaining 550 of the original 986 daily mea- 
surements, each made with 5 different window widths. 
Finally, we corrected the diameter values for several 
geometrical sources of systematic error, and projected 
all measurements to a standard seeing width of 6 fl , 
which is the most frequently observed value. 

Figure 1 shows the unweighted average of these 
diameter values, reduced to a Sun-Earth distance of 
1 AU and plotted against the FFTD window width. 




10 20 30 40 50 

Window width (pixels) 



Fig. 1. — Observed and inferred diameters (i.e., sep- 
arations between limb positions), as a function of win- 
dow width. Observed values are shown as filled cir- 
cles, calculated values for Model 1 as stars and Model 
2 as diamonds. Seeing with a FWHM of 6 Tt was as- 
sumed. Also shown are the results of linear least- 
squares fits to the limb positions, defining the extrap- 
olation to zero window width. 

Also shown are the results of applying the FFTD to 
the seeing-blurred limb-darkening functions derived 
from two different model solar atmospheres (see be- 
low). The agreement between theory and observation 
is satisfactory, although residual differences affect the 
necessary extrapolation to zero a, and are a signifi- 
cant source of systematic error. 

Figure 2 shows the measured time series of diame- 
ter measurements D (projected to a = by linear 
extrapolation) and of a, averaged over Carrington 
rotation periods of 27.275 days. The error bars are 
standard deviations of the mean for each rotation, 
estimated from the dispersion of the daily measure- 
ments within that rotation. The scatter among the 
daily diameter measurements is about 0.4^ rms, and 
arises mostly from time-dependent motions of the so- 
lar image related to atmospheric seeing. The diame- 
ter was essentially constant throughout the 6.3-year 
observing period, aside from a possible but poorly- 
sampled upturn of about 0.1 1i during 1987. The limb- 
darkening function slope, however, varied with a time 
scale of a year or more over a total range of about 
2.5%, being steeper at solar minimum than at maxi- 
mum. 

We estimate Do as the unweighted mean of the 
550 daily values, with a random error equal to the 
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standard deviation of the mean. These values are 

D = 1919.359* ±0.018* . (2) 
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Fig. 2. — (a) Diameters extrapolated to zero window 
width, averaged over intervals of one Carrington ro- 
tation, shown as a function of calendar date, (b) The 
limb-slope parameter a, averaged in the same fash- 
ion as in (a). The solid and dashed horizontal lines 
indicate the slopes found for Models 1 and 2. 

3. Modeling 

As mentioned in the Introduction, the quantity to 
be measured is the radius of the surface at which the 
local temperature equals the solar effective tempera- 
ture. Evidently, this radius is related to Do in Eq. 
(2), but the relationship is not a simple one; it de- 
pends upon the radiation transfer in the outer solar 
atmosphere, and upon the behavior of the FFTD limb 
definition. To infer the correct radius from the obser- 
vations, one must use a physically-based model of the 
solar atmosphere to calculate the emergent intensity 
as a function of distance from the center of the solar 



disk, and then compute the location on this bright- 
ness profile that would be identified as the edge by the 
FFTD. We calculated the limb intensity by integrat- 
ing the equation of transfer along rays through an 
assumed spherically symmetrical solar atmosphere. 
Since the observations were carried out in a relatively 
narrow wavelength region around 800 nm, we consid- 
ered simply the monochromatic continuum intensity 
at this wavelength. 

We have used two models of the solar atmosphere. 
Model 1 was kindly computed by Rodney Medupc 
with the ATLAS9 code (Kurucz 1993). Model 2 was 
obtained as an average of a hydrodynamical simu- 
lation of convection in the upper part of the solar 
convection zone and lower atmosphere, as described 
by, e.g., Stein & Nordlund (1989) but with updated 
physics (Trampedach 1997); the average was per- 
formed at constant monochromatic optical depth, at 
800 nm. The opacity was computed from the ATLAS 
data in both cases. For Model 1 the source function 
S\ was obtained from the ATLAS code, and hence al- 
lowed for mild departure from LTE. For Model 2 LTE 
was assumed, so that S\ = B\, the Planck function. 

To simulate effects of seeing, we convolved the in- 
tensity with a gaussian, with full width at half max- 
imum (FWHM) specified in arc sec and converted to 
linear distance at 1 AU. After this convolution, we 
integrated the intensities over pixels corresponding to 
1 " at 1 AU, to match the observed intensities. We 
folded the pixel-weighted intensities with the FFTD 
weights over the five different windows described in 
section 2, and carried out the subsequent analysis to 
determine the limb position, through extrapolation to 
zero window width, as for the observations. 

The results of applying the observational proce- 
dure to the computed pixel-averaged intensities for 
Models 1 and 2, assuming 6 1 * seeing are shown in Fig. 
1. The observed variation of diameter with window 
width and that calculated from the models agree well, 
especially for Model 1. Of particular interest are the 
extrapolations to zero window width, corresponding 
to the observed limb position measured relative to 
the nominal photosphere of the models; we obtain 
0.47950 Mm and 0.51634 Mm for Models 1 and 2. 

We have tested the sensitivity of the results to var- 
ious assumptions in the calculation. Replacing the 
true source function S\ by B\ for Model 1 changes 
the limb position by much less than 0.001 Mm; thus 
the assumption of LTE in Model 2 is not a significant 
source of error. Using an assumed seeing of less than 
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6 1 * changes the limb position by less than 0.01 Mm, 
confirming the insensitivity of the FFTD to effects 
of seeing. However, the difference between the two 
models obviously remains a source of some concern. 

4. Results and discussion 

We adopted the modified IAU (1976) value of 
1.4959787066 xlO 5 Mm (Astronomical Almanac, 1997) 
for the astronomical unit, and adjusted this value 
by -4.678 Mm to account for the mean displacement 
between the telescope's noontime location and the 
Earth's center, and by +0.449 Mm for the displace- 
ment of the Sun's center relative to the barycenter of 
the Earth-Sun system. This distance, combined with 
D from Eq. (2), yields the Sun's apparent radius. 
Applying the model corrections described in the last 
section, we obtain 

R Q = (695.5260 ± 0.0065) Mm for Model 1 

R Q = (695.4892 ± 0.0065) Mm for Model 2 

We estimate the modeling errors to be l/y/2 of 
the difference between these estimates, or about 0.020 
Mm. Based on the uncertainties in the geometric cor- 
rections that were made to the measured radius, we 
estimate the systematic errors in the measured value 
to be 0.015 Mm, or about twice as large as the ran- 
dom errors. Averaging our results for Models 1 and 
2, and adding the various error sources in quadrature, 
we arrive at our final estimate of 

R Q = (695.508 ± 0.026) Mm 

The inferred solar photosphcric radius is smaller 
by about 0.5 Mm than the normally used value of 
695.99 Mm (Allen 1976). A review of recent obser- 
vations was given by Schou et al. (1997), concluding 
that these were consistent with an angular diameter 
of 1919. 26 ? ± 0.2 1 *, corresponding to Allen's value of 
Rq. This is also consistent with the observed value 
obtained here (cf. eq. 2). However, it appears that 
the observations considered by Schou et al. refer to 
the inflection point of intensity (or, in one case, to an 
FFTD determination) and hence do not contain the 
correction to photospheric radius. Such a correction, 
taking into account the observational characteristics, 
is an essential part of the radius determination. 

Some confirmation of the reliability of the model- 
ing comes from the comparison in Fig. 2 of computed 
and observed slopes of the limb position as function 



of the scan widths. Nevertheless, it is striking that, 
as indicated by the difference between Models 1 and 
2, the major uncertainty in R & appears to come from 
the modeling. Indeed, it is evident that the real solar 
atmosphere is substantially more complicated than 
the one-dimensional model resulting from the ATLAS 
code or the mean model obtained from the hydro- 
dynamical simulations. A more accurate determina- 
tion of the radius correction can probably be obtained 
from a detailed calculation of the limb intensity, tak- 
ing into account the inhomogeneous nature of the rel- 
evant layers, on the basis of the simulations. Such 
an investigation is beyond the scope of the present 
paper, however. 

We find no significant variation in the observed di- 
ameter during the observation period (cf. Fig. 2); an- 
nual averages of the radius for the years 1981 to 1987 
all agree within their measurement errors of ±.037 
Mm. These limits are substantially smaller than di- 
ameter changes reported previously for the same in- 
terval of time (e.g. Ulrich & Bertello 1995, Laclare et 
al. 1996), but arc in agreement with measurements 
by Wittman (1997). On the other hand, the limb- 
position slope shows fairly substantial variations. We 
also note that during solar maximum, the daily slope 
values tended to be highly variable as well as small in 
magnitude; this suggests that the long-term variation 
may result from localized activity-dependent features 
such as faculae. It is plausible that the previously in- 
ferred variations in solar diameter with solar activity 
is in fact a reflection of such variations in the limb- 
darkening slope. 

It is interesting that the value of Rq obtained here 
is somewhat smaller than that inferred from the so- 
lar f-mode frequencies, indicating additional contri- 
butions to the differences between the observed and 
model values of these frequencies. This issue, and the 
effects of the reduction of the model radius on the 
helioseismically determined structure of the solar in- 
terior will be considered elsewhere. We note, however, 
that Antia (1998) and Schou et al. (1997) found sig- 
nificant effects on the helioseismically inferred sound 
speed from corresponding radius changes. 
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